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Abstract 

We investigate the inactive-active phase transition in an array of additive (exclusive- 
or) cellular automata under noise. The model is closely related with the Domany- 
Kinzel probabilistic cellular automaton, for which there are rigorous as well as nu¬ 
merical estimates on the transition probabilities. Here we characterize the critical be¬ 
havior of the noisy additive cellular automaton by mean field analysis and finite-size 
scaling and show that its phase transition belongs to the directed percolation univer¬ 
sality class of critical behavior. As a by-product of our analysis, we argue that the 
critical behavior of the noisy elementary CA 90 and 102 (in Wolfram’s enumeration 
scheme) must be the same. We also perform an empirical investigation of the mean 
field equafions fo assess their quality and find thaf away from the critical point (but 
not necessarily very far away) the mean field approximations provide a reasonably 
good description of fhe dynamics of fhe PCA. 
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1 Introduction 


Cellular automata (CA) are discrete-space, discrete-time synchronous (states change simulta¬ 
neously everywhere) deterministic dynamical systems that map symbols from a finite alphabet 
into the same alphabet according to local (finite range) rules. The origins of CA date back at 
least to the early days of the modem computer era (~1940-1950), when they were conceived as 
model systems for simple self-reproducing, self-repairing organisms and, by extension, logical 
elements operating in parallel and memory storage devices, among others [1, 2, 3, 4, 5] 

Probabilistic cellular automata (PCA), in turn, are CA that evolve by local rules that de¬ 
pend on the realization of some random variable [6, 7, 8, 9] In computer science, PCA are 
used, e.g., to model the impacts of faulty elements (like reading heads, logical circuitry, or neu¬ 
rons) on the computational capability of CA. Besides serving as model systems for the analysis 
of computation—both applied and theoretical, digital or biological—under noise, PCA have 
also been playing a significant role in the elucidation of some deep issues in equilibrium and 
nonequilibrium statistical mechanics [10, 11, 12, 13, 14] We also mention that, currently, CA 
and PCA are very useful tools in the modeling of biological and ecological complex systems 
[15, 16, 17, 18] 

In this work we investigate the very simple noisy additive (exclusive-or) PCA with respect 
to its inactive-active phase transition. The rigorous determination of the parameters for which 
the noisy additive PCA is ergodic is an open problem [19] In the jargon of statistical mechanics, 
this question amounts to determining whether the noisy additive PCA displays a phase transi¬ 
tion between a phase with a single stationary state devoid of active cells (the inactive phase) 
and a phase where the stationary state has a finite density of active cells (the active phase) and, 
if yes, to estimate the critical parameters of the model. We partially answer this question by 
showing that the noisy additive PCA corresponds to the Domany-Kinzel PCA [10, 1 1] in one of 
its parameter subspaces. We have organized this paper as follows. In Sec. 2, the noisy additive 
and the Domany-Kinzel PCA are described, and we show how they are equivalent to each other 
and to other known CA and PCA in the literature. In Sec. 3 we provide simple analytical ap¬ 
proximations to the critical point of the noisy additive PCA by standard mean field techniques 
and in Sec. 4 we supplemenf fhe analysis wifh Monte Carlo simulafions and finife-size scaling 
analysis, obtaining relatively precise estimates for its critical point and critical exponents. In 
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Sec. 5 we make an attempt to assess how good mean field approximations are for the noisy 
additive PCA by measuring (numerically) the discrepancies between some exact quantities and 
the same quantities in the mean field approximation. Finally, in Sec. 6 we summarize our 
resulfs and suggesf a few direcfions for furfher investigation. 

2 The noisy additive cellular automata 

2.1 The p-XOR PCA 

The noisy addifive cellular aufomafon is a fwo-sfafe PCA wifh sfafe space given by Q.\ = 
{0,1}^, wifh A C Z a finife array of |A| = L ^ 1 cells under periodic boundary conditions 
i + L = i, and fransifion funcfion <I>: thaf given fhe sfafe x’ = (xj ... ,x^) of fhe 

PCA af insfanf t G N defermines fhe sfafe = 0(x;,x-^j) of fhe PCA af insfanf 

t + \ according fo fhe following rule: wifh probabilify p G [0,1], xj^^ =A® A+\ exclusive- 
or operation), ofherwise x^^' = 0. In terms of conditional probabilities W (x,-1 x,-,x,+i), we have 
the set of transition rules IT(1|0,0)=0, IT(l|0,l) = IT(l|l,0)=p, and IT (111,1) = 0, with 
lT(0|-,-) = l—IT(1|-,-). This model was first considered by Vasershtein back in 1969 [7], 
among other models that were later to be reintroduced in the literature. We have also considered 
this PCA before, where the phase transition and critical parameters of a closely related asyn¬ 
chronous interacting particle system were established by a transfer matrix-like technique [20]. 
We dub the noisy additive PCA the p-XOR PCA. Its rule table over a symmetric neighborhood 
of radius 1 appears in table 1 . 

The p-XOR PCA interpolates between the trivial CA = 0 at p = 0 and the 

elementary CA 102 (in Wolfram’s enumeration scheme [2]) at p = 1, being a mixture of the 
two (a dilution” of CA 102) for 0 < p < 1. Such probabilistic mixtures of deterministic CA 
have been considered before [21, 22, 23]; in the notation of [23], the p-XOR PCA becomes 
the pl02-(l — p)0 PCA. If p is small, we can reasonably expect that the p-XOR PCA will 
eventually converge to the absorbing state 0 = (0,0,... ,0) devoid of active cells, the corre¬ 
sponding invariant measure being denoted by (5o. It is not known, rigorously, whether there 
exists a critical value p* such that for p > p* the invariant measures of the p-XOR PCA be¬ 
come translation-invariant convex combinations of the form a^o + (1 ~ with 0 < a < 1 
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Table 1: Rule table for the noisy additive PCA. The first row lists the initial neighborhood and the other 
two rows give the probability at which the central cell reaches the state given in the leftmost column. 


111 

110 

101 

100 

oil 

010 001 000 

0: 1 

1-p 

1-p 

1 

1 

1-p 1-p 1 

1: 0 

P 

P 

0 

0 

p p 0 


and fjip the measure that puts mass on configurations with density 0 < /ip(l) < 1. 

We will show that the p-XOR PCA is closely related with the Domany-Kinzel (DK) PCA 
[10, 11], which displays several inactive-active-type phase transitions and for which there are 
rigorous as well as numerical estimates on the transition probabilities. In the next section we 
describe the DK PCA and how the two PCA are related to each other and to other known CA 
and PCA in the literature. 

2.2 The Domany-Kinzel PCA 

The Domany-Kinzel (DK) PCA is a two- (sometimes three-) parameter PCA originally intro¬ 
duced to investigate the relationship between r/-dimensional nonequilibrium models and d + \- 
dimensional equilibrium models, casting some light on this difficult subject [10, 11, 12, 13, 14, 
19]. Its rule table appears in table 2. The two-dimensional parameter space of the DK PCA 
encompasses several well-known, archetypal models in discrete mathematics and theoretical 
physics, among them the directed site {p\ = p 2 ) and the directed bond (p 2 = 2pi — p\) perco¬ 
lation processes, the exactly solvable compact directed percolation process, which is equivalent 
to the zero temperature Glauber-Ising model, the voter model, and diffusing-annihilating ran¬ 
dom walks {p\ = 1/2, p 2 = 1), and a couple of other models, like the elementary CA 90 
(pi = I, P 2 = 0)^ and Stavskaya’s PCA {p\ = p 2 = I — £, the noise parameter of Stavskaya’s 
PCA) [10, 11, 12, 13, 14, 19, 24, 25, 26, 27]. 

From tables 1 and 2, we see that the p-XOR and the DK PCA do not immediately relate 
to each other. Rule tables 1 and 2 are compatible only at the point p = p\ = P 2 = 0, that 
corresponds to the trivial CA 0. Indeed, in the notation of [23], we have that the p-XOR 

^Refs. [26, 27] incorrectly state that the line p 2 = 0 of the DK PCA corresponds to the diluted elementary CA 18. 
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Table 2: Rule table for the Domany-Kinzel PC A. Same notation as in table 1. 


111 

110 

101 

100 

Oil 

010 

001 

000 

0: l-p2 

1 -pi 

1 -P2 

1 -pi 

1 -pi 

1 

1 -Pl 

1 

1: p2 

Pi 

Pi 

Pi 

Pi 

0 

Pi 

0 


PCA is the /tl02-(l — p)0 PCA, while the DK{pi,p 2 ) PCA becomes, at special points, PCA 
DK(p,0) = p90-(l — p)0, DK{p,p) = p250-(l — p)0, and DK(0,p) = pl60-(l — p)0, and 
each of these PCA does not relate with PCA pl02-(l — p)0 except at the trivial point p = 0. 

Note, however, that the dynamics of the DK PCA occurs in two separated sublattices 
that do not interact, each of which evolves by a rule that can be written as 
see figure 1. If we recast the DK PCA in one of its sublattices as a PCA over a symmetric 
neighborhood of radius 1 again, we obtain the rule table given in table 3. Now, comparing 
the DK PCA dynamics in one of its sublattices with the p-XOR PCA dynamics, we see that 
the p-xor PCA corresponds to the DK PCA on the line p 2 = 0. On this line, it is known that 
the “diluted” CA 90 and the DK PCA coincide and that they display an inactive-active phase 
transition at p* ~ 0.81 [10, 11, 28, 29, 30, 31, 32]. 

We could have simplified our discussion above by looking af fhe rules of fhe p-XOR and 
fhe sublaffice DK PCA in fheir asymmefric versions and noficing fhaf fhey 

are equal when p 2 = 0. Buf fhe discussion using symmefric neighborhoods of radius 1 indicates 
fhaf fhe critical behavior of fhe diluted CA 90 (or PCA p90-(l — p)0), represented by fhe DK 
PCA af fhe poinf p 2 = 0, and fhe diluted CA 102 (or PCA pl02-(l — p)0), represented by 
the p-XOR PCA, must be the same. The fact that these two CA/PCA may display the same 
critical behavior seems to have gone unnoticed in the literature so far. Note, however, that it has 
been found that some properties may differ whether one uses the full lattice or the sublattice 
version of the DK PCA [33]. It should be remarked, however, that for the cases investigated in 
[33] there is an accompanying change in the neighborhood that changes the nature of the rules 
(the full lattice case includes the cell itself in the rule); this can be clearly seen if we blow the 
resulting rules up to a symmetric neighborhood of radius 1, as we did for the p-XOR PCA—we 
would then discover that the two models investigated there are compatible only on the trivial 
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Figure 1: The dynamics of the DK PC A occurs in two separated sublattices that do not interact. In the 
figure, the circles together with the dashed arrows display one of the sublattices. 

Table 3: Rule table for the Domany-Kinzel PCA in one of its sublattices recast as a PCA over a sym¬ 
metric neighborhood of radius 1. Same notation as in tables 1-2. 


111 

110 

101 

100 

Oil 

010 

001 

000 

0: l-p2 

1 -pi 

1 -pi 

1 

1 -P2 

1 -Pl 

1 -Pl 

1 

1: p2 

Pi 

Pi 

0 

P2 

Pi 

Pi 

0 


point pi= p 2 = 0 . 


3 Mean field analysis 

The Markovian dynamics of the probability distribution Pt{x) of the states x G Q.j^ of the PCA 
is governed by the equation 

P,+i{x')='£Wix'\x)Prix), ( 1 ) 

X 

where the summation runs over all x G and IT (at' | a;) ^ 0 is the conditional probability 
for the transition at —)■ at' to occur in one time step. When the cells of the PCA are updated 
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simultaneously and independently we have 


L 

W{x'\x) = Y\yVi{x'i\x), with 'Y_yVi{x\\x) = \. 


(2) 


For the p-XOR PCA, Wi{x\ \ x) = W{x\\xi,Xi+\), independent of i —the PCA is homoge¬ 
neous in space. The marginal probability distribution Pt{xi,... ,x„) of observing a block of n 
consecutive cells in state (.xi,... ,x„) is obtained from Pt{x) by summing it over the variables 
Xn+i, ■■■,xl- From (1) and (2), we derive that the dynamics of Pt{xi,... ,Xn) is given by 




^ 1 , 




Pt{xi,...,x„+i). 


(3) 


We see that the probability of observing n consecutive cells in a given state at instant f + 1 
depends on the probabilities of observing the state of n + 1 cells at instant t. 

To proceed with the calculations, either we solve the full set of L coupled equations (3) 
exactly or truncate the hierarchy of equations (3) at some point to get a closed set of equations. 
The simplest approximation is obtained by taking 


Pt{x\,X2) «Pf(xi)Pf(x2), 


(4) 


while higher order approximations can be obtained by the generalized approximation 


Pt{x\,...,Xn) 


Pf(xi,...,X„-i)Pf(x2,...,X„) 

Pf(x2,...,X„_i) 


(5) 


This approximation has been described in the context of cellular automata in [34]. For prag¬ 
matic expositions of the technique see [35, 36, 37, 38, 39, 40, 41], while further developments 
and applications can be found in [39, 41, 42, 43]. 

We note the following features of the p-XOR PCA and equations (3) and (5): 


(/) Since W{ - \x\,X 2 ) = W(- \x 2 ,xi), the PCA is reflection-symmetric and the marginal 
probability distributions observe P(xi,X 2 ,... ,x„) = P(x„,x„_i,... ,xi); 


(//) The cluster approximation (5) is consistent with the reflection-symmetry of the marginal 
probability distributions, i. e., the approximations of P(xi,X 2 ,... ,x„) andP(x„,x„_i,... ,xi) 
according to (5) coincide; 
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(Hi) If we take P(l, 1,..., 1) on the left-hand side of (3), only the terms P(1,0,1,...) and 
P(0,1,0,...) appear on the right-hand side, because both 1^(110,0) and Ik(111,1) are 
zero. Moreover, for even number or arguments P(1,0,..., 1,0) = P(0,1,... ,0,1) by 
property (/) above, further simplifying the mean field equations. 

We now proceed with the mean field approximafion fo fhe full probabilify disfribufion of con¬ 
figuration of fhe p-XOR PCA up fo fhe level of fwo cells. 


Single-cell approximation. From fable 1 and equafions (3)-(5), fhe single-cell approxima¬ 
tion for P+i (x = 1) reads 


P,+i(l)=pP,(0,l)+pP,(l,0) = 2pP,(0,l) 
«i2pP,(0)P,(l) = 2p[l-P,(l)]P,(l). 

In fhe sfafionary sfafe, Pt+i{ 1) = Pf (1) and we obfain 


( 6 ) 


P(l)=2p[l-P(l)]P(l), 


(7) 


wifh solutions P(l) = 0 and P(l) = {2p — \)/lp. The firsf solution corresponds fo fhe single 
sfafe 0 = (0,... ,0) devoid of acfive cells. The ofher solution is negafive (unphysical) for p < 
1 /2 and posifive for p > 1 /2, fhus predicfing a phase fransifion af p* = 1 /2 from fhe inactive 
sfafe fo a sfafe wifh densify 0 ^ P(l) = (2p — l)/2p ^ 1/2 of active cells. Note fhaf in fhis 
approximation, P(l) = 1/2 af p = 1, which coincides wifh fhe exacf value of fhe sfafionary 
densify of acfive cells of CA 90 [2], which is fhe DK PCA af pi = 1, p 2 = 0. The approximation 
P(xi,X 2 ,X 3 ) Ri P(x\)P{x 2 )P{x-i) of [2, (3.4)] also estimates fhe sfafionary value P(l) = 1/2 for 
CA 102 (recall fhaf fhe p-XOR PCA is, in fhe nofafion of [23], fhe mixed PCA pl02-(l —p)0). 


Two-cell approximation. Since P,(0,1) -|-P,(1,1) = P((l) and Pr(l,0) = Pf(0,1), we only 
need to write equations for Pf (1) and Pj (1,1) to characterize the two-cell mean field approxi¬ 
mation of the p-XOR PCA. The equation for P^ (1) is (6), and from table 1 and equation (3) we 
obtain 

P,+i(l,l)=p2p,(0,l,0)+p2p,(1,0,1). (8) 
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The approximation (5) provides 


P, (0,1,0) 


P,(1,0,1) 


P,(0,1)P,(1,0) 

^r(l) 

P,(1,0)P,(0,1) 

m 


such that, in the two-cell approximation, 


^m(l,l) 


2 [P,(1)-P.(1,1)]^ 

P,(l)[l-P,(l)] 


(9) 


( 10 ) 


We want to solve for P(l), the order parameter” of the model. From the exact part of equation 
(6) we must have P(l, 1) = {2p — l)P{l)/2p in the stationary state. Imposing the stationarity 
condition /’m(l,l) = P,(l, 1) in (10) and plugging in the expression for P(l, 1) in terms of 
P( 1) furnishes the solution 

= ( 11 ) 

for the density of active cells in the stationary state of the p-XOR PCA in the two-cell mean 
field approximation. This result predicts an inactive-active phase transition at the critical point 
p* = 2/3; again, we have P(l) = 1/2 at p = 1. 

The two-cell mean field approximafions for the dynamics of the p-XOR PCA provide the 
lower-bound p* ^ 2/3 for the critical value of the model. Upper bounds are more difficult 
to obtain [24]. Although it is possible to improve the lower bound by pushing the mean field 
approximation to higher orders, the equations soon become somewhat ungainly. We shall pro¬ 
ceed to Monte Carlo simulations and finite-size scaling analysis to gain a better understanding 
of the phase transition of the model. 


4 Finite-size scaling analysis 

Our Monte Carlo simulations of the p-XOR PCA ran as follows. For a given p, the PCA is 
initialized with each cell x/ = 1, 1 ^ ^ L, drawn independently with probability 1/3, i.e., 
~ Binomial(L, 1 /3). Stationary state quantities, e.g. the density of active cells Pl = P{1) = 
are then sampled after the system is relaxed through 5L Monte Carlo steps (MCS), 
with one MCS equivalent to a synchronous update of the states of all L cells of the automaton. 
This amount of relaxation proved sufficient to reach the stationary state away from the critical 







Figure 2: Stationary density P(l) of active cells in the one- and two-cell mean-field approximations 
together with Monte Carlo simulation data for a lattice of L = 16000 cells initialized randomly with 
density ~l/3 (cf. Sec. 4). Each symbol in the Monte Carlo curve is an average over 2500 samples. 
These data estimate the critical point at p* ~ 0.81. 

point. Note, however, that except for the data in Fig. 2, our results were obtained from time- 
dependent simulations, not from averages in the stationary state. 

According to the general theory of critical phenomena for equilibrium and nonequilibrium 
systems [25, 26, 44], we can assume that close to the critical point the density pi of active 
cells obeys the finite-size scaling relation 

PL(f;AL) ~t-^/''ll<I>(ALt'/v|i^fVx/v||/^)^ (12) 

with = p — 0 and L the size of the PCA array. For a very large system, pi{f,Ai) ~ 

<I)(Ai,f^/''il), with <I>(v <C 1) ~ const and <I>(x 1) ~ The investigation of the time- 

dependent profiles pL{t',Ai) fhen allows fhe simultaneous determination of p\ and 8 = j8/V||. 
Perusal of (12) and derived relations furnish the other exponents, see, e.g., [26]. 

Figure 2 displays the density profile Pi for an aufomafon of L = 16000 cells in fhe sfafion- 
ary slate. The steep Iransilion aboul p£ ~ 0.81 anticipates a small value for fhe exponenl j8. 
Note fhaf fhis rough estimate of p* agrees well wifh fhe critical poinf 0.80 < p] < 0.81 found 
before for fhe DK PCA phase Iransilion along fhe line P 2 = 0 [10, 11, 28, 29, 31, 32, 41]. 

To estimate p* more precisely, we plof Ptit) close lo p ~ 0.81 for some large L. On fhe 
critical poinf, ptit) ~ and we can estimate 5 by plotting log^[pL(t/f?)/pz,(t)] againsl 1/t 
for some small b. Our dafa for L = 16000 and b = \0 appear in figure 3, where each curve 
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Figure 3: (a) Logarithmic plot of PL(t) with L — 16000 cells, (b) Instantaneous values of d obtained 
from the curves in (a). In panel (a) we have, from the lowermost curve upwards, p = 0.80930, 0.80934, 
0.80938, 0.80942, 0.80946, and 0.80950, while in panel (b) this order is reversed. The dashed line in 
panel (b) indicates the best value available for 5 dp = 0.159464(6). 

is an average over 2500 realizations of the process. From these data we extract the estimates 
p* = 0.80938(4) and 5 = 0.159(2), where the numbers in parentheses indicate the uncertainty 
in the last digit of the data. The value for 5 was obtained from extrapolations of the curve at p = 
0.80938, which give 5{t) = —12.746t^^ +0.1593 with a correlation coefficient = 0.1506, 
or, alternatively, d{t) = 0.1593 exp(—80.5t^^) with a correlation coefficient = 0.1516. 

The exponent V|| can be obtained by plotting pi{t) against and tuning V|| to achieve 
data collapse with different A^. Figure 5 shows the collapsed curves obtained with p* = 
0.80939, 5 = 0.159, and V|| = 1.75. We could not discern the value of V|| more precisely 
than by ±0.05. Otherwise, the data collapse is very sensitive to p* and we used this fact to fur- 
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Figure 4: Data collapse of the scaled time-dependent density profiles on the critical point p* = 0.80939 
for some 4000 P L P 16000. Best collapse was obtained with 5 = 0.158 and z = 1.55. 


ther bracket the critical point to p* = 0.80939(3). Combining 5 = 0.159(2) and V|| = 1.75(5) 
furnishes j8 = 5v|| = 0.278(9). 

The third independent exponent can be obtained by plotting versus t/U- for differ¬ 

ent L and tuning z until data collapse for some z. Since z = V|| /by definition, this procedure 
also gives Vj^, once V|| is known. The finite-size curves appear in figure 4. We found besf dafa 
collapse wifh z = 1.55(5) and 8 = 0.158. 

Exponenfs 5, V||, and z suffice fo defermine fhe universalify class of critical behavior of 
fhe model, fhe ofher exponenfs following from hyperscaling relations [25, 26]. The besf values 
available for fhe crifical exponenfs of fhe (1 + 1)-dimensional direcfed percolafion process on 
fhe square laffice are 5 dp = 0.159464(6), ypp = 1.733 847(6), jSdp = 0.276486(8), andzop = 
1.580745(10) [45, 46]. Thus, wifhin fhe error bars, our estimates of fhe exponenfs of fhe p- 
XOR PCA—namely, 5 = 0.159(2), V|| = 1.75(5), j3 = 5v|| = 0.278(9), and z = 1.55(5)—pul 
ils phase Iransilion in fhe (1 + 1)-dimensional directed percolafion universalify class of crifical 
behavior. 


5 An assessment of the mean field approximation 

While in fhe one hand if is hard fo obfain rigorous bounds on fhe dislance befween fhe full exacf 
invarianf measure for fhe PCA and fhe approximalions provided by fhe mean field equafions. 
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Figure 5; Data collapse of the scaled time-dependent density profiles for several different values of p. 
The lower (upper) branches correspond to p <p* (p> p*). The best data collapse was obtained with 
p* = 0.80939, 5 = 0.159, and V|| = 1.75. 


on the other hand one can worry that the one- and two-cell mean field approximations are not 
very good away from p = 1, as it is clear from figure 2. Moreover, alfhough we can reasonably 
expecf fhaf higher order mean field approximafions become heller and better, if is our experi¬ 
ence lhal Ihe convergence of Ihe approximations lo Ihe limil values close lo Ihe critical regions 
is slow wifh fhe order of Ihe approximalion, making simple exfrapolalion slralegies nol very 
useful—Ihe facl fhe slafionary measure becomes singular (in fhe limil L —)• oo) in the vicinity 
of a phase transition only makes things tougher. 

In this section we take an empirical look at the ratio between the k-cell mean field approx¬ 
imation {xi,...,Xn) and fhe exacl marginal probabilily dislribufion Pi{xi,... ,Xn), 

Plllixu-.-^Xn) 


Rf\xi,...,Xn) = 


Pf+l(xi,...,X„)’ 


(13) 


as a proxy fo fhe “qualify” of fhe mean field approximation. The index t + 1 on fhe righl-hand 
side of (13) is a matter of convenience; we can look af r\’^\x\,. .. ,x„) as a predictor for fhe 
discrepancy of fhe probabililies being measured in fhe nexl step of fhe dynamics of fhe PCA. 
Clearly, R^^ (xi,... ,x„) = 1 for k > n, since in Ihis case p/^^ (xi,... ,x„) = P,(xi,... ,x„). Note 
lhal R^f^\x\, ...,Xn) is related wifh fhe n-poinl equal-lime connected correlation function [44] 
gI‘^\xi,. .. ,Xn) = [Pt{x\) — P{x\)\ ■ ■ ■ [Pt{xn) — P(x„)], bul wc will nol explore Ibis relation here. 
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Figure 6: Ratios (upper panel) and Rp^(l,l) (mid panel) for several values of p. Time t is 

given in Monte Carlo steps. In the upper panel we have, from the uppermost curve down, p = 0.77 and 
0.79 {p < p*, dashed lines), 0.81, 0.83, 0.85, 0.87, and 0.89 (p > p*, solid lines). In the mid panel this 
order is reversed. The lower panel displays the stationary values reached by R^pl) and Rppl, 1) for 
different values of p. 
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From equations (6), (8) and (9) we obtain 




mm) 

Pt{0,l) 


(14) 


for the discrepancy of the probability density P? (1) in the single-cell approximation and 

P(0,1)P(1,0) 

' ^ ^ P(0)P(1)[P(0,1,0)+P(1,0,1)] ^ ^ 

for the discrepancy of the probability density Pf(l, 1) in the two-cell approximation. Note that 
Pp^(l) and Pp^(l, 1) are defined only for p ^ p*, because for p < p* we have P(l) = 0, and 
so every other marginal P(..., 1,...) = 0 as well. We define p(^)(l) = 0 and P^^^(l, 1) = 0 for 

p^p*. 

To compute (14) and (15), we measured the probabilities Pf(l) = 1 — P/(0), Pi(0,1) = 
P,(1,0), Pr(0,1,0), and P/(1,0,1) by Monte Carlo simulations. For each value of p, quantities 
were averaged over 2500 realizations of the dynamics. We did not find significant differences 
in the data for different sizes of the array (L = 2000, 4000, and 8000 cells), so we report results 
for L = 4000 only. Our results appear in figure 6. 

At t = 0, Pq^^( 1) = 1 and Po^^(l, 1) = 1, because the initial state is an uncorrelated 
random state with Po(xi,... ,x„) = Po(xi) • • •Po(x„). Otherwise, we observe rapid convergence 
of Pp^(l) and Pp^(l, 1) to their stationary values when p is away from p*\ close to p* correla¬ 
tion lengths increase and it takes longer for these quantities to reach their stationary values. We 
observed the empirical bounds p|^^(l) <1.4 and Pp^(l, 1)^1 for all p; if we consider only 
the region p > p*, then the empirical bounds read 1 ^ Pp^(l) <1.4 and 0.6 < Pp^(l, 1) ^ 1. 
It is thus clear that p/'^(l) is an upper bound to P((l), while P/^^(l, 1) is a lower bound to 

Our conclusion is that away from the critical point (but not necessarily very far away), 
the mean field approximations provide a reasonably good description of the dynamics of the 
PCA. In particular, even the relatively crude one-cell approximation for the density of active 
cells is never far off the value of the empirical density by more than ~ 40%, while the two-cell 
approximation falls even closer to the empirical value of its respective marginal probability 
density. For example, the determination of the critical point from the two-cell approximation 
differs from the empirical value by ~ (2/3)/0.81 ~ 5/6, i.e., by just ~ 1/6 of the empirical 
value. The short range of the interaction in the p-XOR PCA (and in PCA in general) may be 
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partially responsible for this good agreement already at the level of the two-cell approximation. 
One cannot, however, rely on the mean field approximations to determine critical exponents 
(they are always classical”) or phase diagrams with more than one parameter, in which case 
the loci of the critical manifolds and multicritical points may differ significantly from the ones 
provided by the mean field approximafions [25, 26, 41, 44]. 

6 Summary and conclusions 

We esfimafed fhe critical poinf of fhe p-XOR PCA af p* = 0.80939(3) and found fhaf fhe 
model belongs fo fhe (1 + 1)-dimensional direcfed percolafion universalify class of crifical 
behavior. This value of p* agrees, wifhin fhe error bounds, wifh some (buf nol all) values found 
before for p\ in fhe DK PCA on ifs line P 2 = 0 [10, 11, 28, 29, 31, 32, 41], buf improves 
former estimates by af leasf one order of maginfude—fhe besf value currenfly available for 
fhis quantify was p\ = 0.8087(5) [31]. Our numerical analyses look approximately 1630 CPU- 
hours (developmenf and fesfs discounted) on Power 755 RISC processors running C/GCC code 
af 3.4 GHz over IBM AIX. 

Our idenfificafion of fhe p-XOR PCA as an inslance of fhe DK PCA helps fo make in- 
formalion exislenf abouf fhe later available in fhe sludy of fhe former. If is known for some 
time lhal fhe diluted CA 90 (i.e., fhe p-XOR PCA) and fhe DK PCA are equivalenf on fhe line 
P 2 = 0 of fhe later [30], buf we have furfher idenfified fhe equivalence of fhe p-XOR wifh fhe 
diluled CA 102, which should fhus have fhe same crifical behavior as fhe DK PCA over fhe line 
P 2 = 0 and fhe diluted CA 90. The relafionship befween fhe DK PCA and two-dimensional 
classical spin models and the detailed knowledge of its phase diagram, that displays an ex¬ 
tended chaotic phase, may help to classify the behavior of the p-XOR PCA and its siblings 
according to Wolfram’s classes 1-4 [47] (see, however, a critique of this program in [48]). 

The analysis of the mean field approximations showed lhal fhe clusler approximafion pro¬ 
vided by fhe general prescriplion (5) is acceplable even for cluslers formed by as few as Iwo or 
Ihree cells. Clearly, fhe approximafion is heller fhe farlher one is away from fhe crifical poinf, 
when correlation lenglhs diverge. A differenl approach lhan lhal presenfed in section 5 would 
be fo compare fhe analylical expressions (7) and (11), and possibly higher order approxima¬ 
tions, wifh fheir empirical counlerparls. The case of P(l) can be apprehended direcfly from 
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figure 2. Moreover, the ratio = P(0)P(1)/P(0,1) could perhaps be used to locate the 

critical point more efficiently (less computation, more precision) than with P(l) alone. We 
intend to pursue these lines of inquire further elsewhere. 

Finally, we believe that a more detailed investigation of probabilistic mixtures of deter¬ 
ministic CA with known properties [21, 22, 23] is worth carrying out further. 
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